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, Models of leptogenesis are constrained by the low reheat temperature at the end of reheating 

(N 

associated with the gravitino bound. However a detailed view of reheating, in which the 
maximum temperature during reheating, Tmax, can be orders of magnitude higher than the reheat 
On . temperature, allows for the production of heavy Majorana neutrinos needed for leptogenesis. But 
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CN ■ then one must also consider the possibility of enhanced gravitino production in such scenarios. 
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\^ \ considered gravitino production in the post-reheating radiation dominated era. We find that the 
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In this article we consider gravitino production during reheating, its dependence on Tmax; and 
its relevance for leptogenesis. Earlier analytical studies of the gravitino abundance have only 



gravitino abundance generated during reheating is comparable to that generated after reheating. 
This lowers the upper bound on the reheat temperature by a factor of 4/3. 
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I. INTRODUCTION 



Gravitinos in supersymmetric theories can have important cosmological consequences. Sta- 
ble gravitinos can overdose the universe while unstable gravitinos can affect the expansion 
rate of the universe during eras prior to their decay. The decay products of unstable grav- 
itinos can also overdose the universe or affect light element abundances generated during 
nucleosynthesis. These cosmological consequences are a function of the gravitino energy 
density, Pq = rriQnQ, where ttIq and tIq are the mass and number density of gravitinos. 

In a non- inflationary universe, tIq ~ and therefore cosmological constraints on the 
energy density of grav,tinos provde bounds on and eqnivalently on the sea.e of super- 
symmetry breaking In an inflationary universe, is also a function of the reheat 

temperature, and so for a fixed rriQ, often taken to be 0(100 GeV — ITeV), cosmological 
constraints on the energy density of gravitinos provide an upper bound on the reheat tem- 
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131] the number density of gravitinos is obtained 



by considering gravitino production in the radiation dominated era following reheating. It 
is presumed that = at the beginning of the radiation dominated era. Gravitinos are 
then produced through thermal scattering and is found to be proportional to the reheat 
temperature, T^eh, which is the temperature of the thermal plasma at the beginning of the 
radiation dominated era when the inflaton field has decayed completely and the energy 
density of the universe is dominated by the inflaton decay products. The cosmological 
constraints on n^, then provide an upper bound on Tj-^h of 10^"^ GeV. 

If, as in Refs. i3i,Q,B,y,0,y,y,Q,i [ill . [12I . I13I , one assumes instantaneous reheating after 

inflation or that T^eh is the maximum temperature during reheating, then the upper bound 
on Treh makes it difficult to create sufficiently high number densities of GUT gauge and Higgs 
bosons whose decays could generate the baryon asymmetry of the universe. Similarly, the 
gravitino bound constrains leptogenesis models in which the lepton asymmetry is generated 
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16|, 



by the decay of heavy bosons or fermions [l^. However, as discussed in Refs. 
after the inflationary era the temperature does not rise instantaneously to T^eh but rises 
initially to a maximum temperature Tmax and then falls to Treh- In Ref. jl^ the authors 
then argue that T^ax can be as high as lO^Treh and that sufficient numbers of the lightest 
heavy Majorana right-handed neutrino of mass ~ lOT^eh can be produced during reheating 
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to allow for successful leptogenesis. This issue has also been studied in Ref. 

While the above scenario considers the possible production of heavy neutrinos during 
reheating it does not consider the possible enhancement in the gravitino production as well. 
If the gravitino abundance generated post-reheating is proportional to the maximum tem- 
perature during that era, namely T^eh, one should ask whether the gravitino abundance 
generated during reheating is proportional to the maximum temperature, Tmax, of the re- 
heating era. If this abundance is larger than the abundance generated in the post-reheating 



era it could affect the viability of the leptogenesis scenario of Ref. [16|. Therefore in this 
article we explicitly calculate the gravitino abundance generated during the reheating era. 
We then compare it to the standard calculation of the gravitino abundance generated after 
reheating in the radiation dominated era, and discuss its relevance for leptogenesis models. 

As argued above, one might expect that the gravitino abundance generated during re- 
heating will be a function of Tmax- Interestingly, what we find is that by manipulating the 
relations between Tmax, ^reh, the infiaton decay rate and the scale of inflation the dependence 
on Tjnax cancels out and the gravitino abundance is proportional to T^eh only. Furthermore, 
while one would not expect the gravitino abundance generated in the reheating and the 
post-reheating eras to be similar that is indeed what we find. The resulting contraint on the 
reheat temperature and on leptogenesis models is then only slightly stronger than before. 

Our results are valid for a reheating scenario that does not include pre heating. Gravitino 



production during preheating has been considered in Ref. |18l . \li 
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II. PRODUCTION OF GRAVITINOS 

During infiation the Universe cools down by several orders of magnitude. Subsequently 
the infiaton decays while performing coherent oscillations about the minimum of its poten- 
tial. Very soon after the infiaton enters the oscillating phase the temperature of the universe 
rises to a maximum value 1^ 

T^ax ^ 0.8g;'/'Mi/' {T^Mp.f , (1) 

where Mi = V^^^\ Vi being the vacuum energy density during the infiationary epoch (taken 
to be constant), g^, is the number of relativistic degrees of freedom and is the decay rate 
of the infiaton field. Subsequently, the temperature of the thermal bath falls approximately 
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as R~^/^ 15|, where R is the scale factor of expansion of the Universe. This particular 
dependence on R goes on until the universe becomes radiation dominated when the inflaton 
ield decays completely at treh = Tl^. The temperature of the universe at treh is given by 
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reh 



a/2 



(2) 



In the following we examine the production of gravitinos during reheating between Tmax and 
Treh, and during the subsequent radiation dominated era, and discuss its consequences. 

Gravitinos are produced by the scattering of the inflaton decay products; see, for example. 
Tables 1 in Refs. 7|, ll2| for a list of processes. The number density of gravitinos generated 
is then given by the solution of the Boltzmann equation 



dt 



+ 3Hn 



G 



{T.tot\v\)n 



(3) 



where n = {({3)/tt'^)T^ (C(3) = 1.20206.. is the Riemann zeta function of 3), Stot is the total 
scattering cross section for gravitino production, v is the relative velocity of the incoming 
particles, and (...) refers to thermal averaging. We have^nored the gravitino decay term 
above as the gravitino lifetime is 10^~®(100 GeV/m(5)s [7] and is not relevant during the 
gravitino production era for gravitinos of mass 10^~^ GeV. We may re-express this equation 
as 

{^tot\v\)n\ (4) 



dT 



G 



In SU{N) supersymmetric models with Uf pairs of fundamental and antifundamental chiral 
supermultiplets, (Stot|'i^|) is given by [3] 



(Stotl^^l) = 



a 

M2 



1 



1 + 



32C(3) 



X [{ln{TymU + 0.3224}(iV + nj) + 0.5781nj] , 



(5) 



where rrig is the gaugino mass and rrig th is the thermal mass of the gauge boson which is 



given as 



:gUN + nf)T' 



(6) 



In the above equations gN{N = 1,2,3) are the gauge coupling constants corresponding 
to U{1)y, SU{2)l and SU{3)c respectively and M = Mpi/Vs^ ~ 2.4 x 10^^ GeV is the 
reduced Planck mass. (For U{1) gauge interactions, N'^ — 1-^1 and N + Uf ^ Uf, where 
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Hf is the sum of the square of the hypercharges of chiral multiplets 29|].) Using the one 
loop /3-function of MSSM, the solution of the renormalization group equation for the gauge 
coupling constants is given by 

-1/2 



9n{T) 



UN_ 

87r2 



HT/Mz) 



(7) 



with 61 = 11, 62 = 1, &3 = —3. It is presumed here that inflaton decays perturbatively 
and the products thermalise quickly as discussed in Appendix A of Ref. jl6|. Also see 



Refs. 
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27( 1 for an alternate description of reheating and the gravitino bound. 



A. Gravitino production during reheating 

If the potential for the oscillating inflaton field is dominated by the mass term, then the 
energy density of the inflaton field scales as during reheating. (We ignore change in 

the inflaton energy density due to decays.) Einstein's equation then implies 

where pmax and -Rmax are the inflaton energy density and scale factor at Tmax- Solving the 
above equation for R we get 



R Rmax 



2/3 

(9) 



where 

3 Mpi 



(10) 



(For t:^t^^^,Rr~. t2/3.) 

For Eq. (jl]) we require T and H as functions of the temperature T. During reheating the 
energy density of the relativistic particles is given by p,. oc i?"^/^ [15|]. Since pr ~ T^, we get 

Using Eq. ([9]) the T — t relation is then obtained as 
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This implies that 



f = --TK 



8 [|-f^max(i — ^max) + l] 



/ T 



~^TH^,^[-^] . (13) 

J- rriF 



The Hubble expansion parameter in terms of T is then given by 



R 8T „ / T ^ ^ 



Using Eqs. ( fT3l) and dH]), Eq. (jl]) can be rewritten as 



- —ri/ 



dT 
where 



-CT, (15) 



Solving for ra^ in the regime T^ax to Treh we get 



C-l^^mX- (16) 



"0(1-,=,,) = f r»u - (17) 

^ \ rch max / 

C 

= — T^eh for Tmax > Treh • (18) 
D 

a has been taken to be constant although there is a log dependence due to the running 
gauge couplings in the one loop (3 function of MSSM. 



B. Gravitino production in the radiation dominated era 

After the infiaton field decays completely at treh the universe enters the radiation dom- 
inated era. Unlike the reheating era during which the entropy continuously increases, in 
the radiation dominated era the total entropy remain constant (except for epochs of out- 
of-equilibrium decays). Therefore it is useful to express the abundance of any species i as 
Yi = n-i/s, where rii is the number density of the species i in a physical volume and s is the 
entropy density given by 

' = '^9*T' . (19) 
= 228.75 in the MSSM. With this definition, Eq. (g]) reads as 

T^ = {j:,M)Yn. (20) 
6 



For the radiation dominated era, 



T = T,eh r , (21) 

[2ifreh(i — ^rch) + 1]^ 



where 



This imphes that T is 



H,., = \r^^. (22) 
\/ 90 Mpi ^ ^ 



" reh 



rp _ "_renrp3 _ / ^*ren" i - 

^~ tTJ ~ V^o~j M- ^^^^ 



Then 

rfy^ _ / 90 1 \/a\/C(3)^^ 



Assuming a to be independent of temperature and integrating the above equation from Treh 
to Tf, the final temperature, we get the number density of gravitinos at Tf to be 



x(^)f^) (25) 



Since most of the gravitinos are generated close to Treh "we have ignored the variation of 
with temperature and used g^ch in the final expression. Using Eqs. ( fTSj) . ( fT6l) and ( fT9|) the 
first term on the right hand side of the above equation is given by 

This term is usually neglected while estimating the gravitino abundance |3|, U, [sl, la, LZl, [sl, Isl, 
' 12I . [1^. However as we see below this is comparable with the second term in Eq. 
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(!25|) . Now using Eq. (126|) in Eq. (125|) we get the effective number density of gravitinos at 
Tf to be 



90 

+ M ( ) T,eh 



4 

max 
9 -ff max Treh 



(27) 



where we have used Tf ^ Treh- Relating T^ax to Treh from Eqs. ([T]) and ([2]) the number 
density of gravitinos in a comoving volume is then given as 

yciTf) = ^ (^-^] ' ( ^—jjA (1-0 + 3.0) , (28) 

where we have used g„eh in the expressions for Tmax- 



III. DISCUSSION AND COMMENTS 



Generalizing Eq. (fT7|) for any T during reheating one sees that does not vary mono- 
tonically during reheating. As shown in Fig. ([T]), rises dramatically from T^ax to 
Ti = and then falls from Ti to Tj-eh- However ~ T^^ and the number density per 

comoving volume, Uq = (l/R^^^) dT d / dT (nQR^) , is proportional to — l/T^^x) 

and so most gravitinos are produced close to T ~ Tj-^h- {s ~ and, as shown in Fig. ([2]), 
Yq increases steadily in the reheating phase for Treh < T < Tmax-) For T < Treh, one can 
show from Eq. fl25|) that dn^/dT is always greater than 0, indicating that the gravitino 
number density is always decreasing during the radiation dominated era. (If Y^(Tj.eh) is set 
to 0, then first increases and then decreases.) Yq is always increasing but for T ^ T^h 
it becomes approximately constant, as seen in Fig. ([2]). 

We now address a concern as to whether one can distinguish production as occuring in 
the reheating era or in the radiation dominated era when production in both eras occurs 
close to Treh- In our analysis, for T < Tmax during reheating the ratio of the radiation 
and inflaton energy densities increases as Prad/Pinf ~ R~^^^/R~^ R^^"^ while the 

gravitino abundance increases as Yq ~ T^^ (from Eq. (12^ ). Therefore much of the change 
in the gravitino abundance in the reheating era occurs when prad ^ Pinf- Similarly, Eq. (!25|) 
implies that in the radiation dominated era the gravitino abundance changes linearly with 
T ~ while the inflaton energy density falls exponentially fast. Again much of the 

gravitino production at T < Treh will occur when it is valid to treat the dynamics of the era 
as due to radiation only. Therefore we believe that our treatment of the problem is valid. 

From Eq. fl25]) it is clear that the gravitino production during the reheating period is 
1/3 of that during the radiation dominated era. A priori one would not have expected the 
gravitino production in both these eras to be similar. Furthermore, the gravitino abundance 
generated during reheating can be re-expressed as independent of Tmax, which is also unex- 
pected. Interestingly, the contribution to Yq from the reheating era is linearly proportional 
to Treh, as it is for the radiation dominated era. It is then straightforward to revise the earlier 
upper bound on Treh from the cosmological constraints on Uq. The bound of 10^~^ GeV will 
now be lowered by a factor of 4/3 and thus is not greatly affected. Since Tn^ax oc "\/Treh, 
Tmax is also not much affected. Therefore heavy particles of mass greater than Treh can 
still be produced during reheating and leptogenesis scenarios are not significantly further 
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FIG. 1: The gravitino number density, Uq, the gravitino number density per comoving volume, Hq, 
and the entropy density, s, are plotted as a function of the temperature T during reheating. Tjch and 
Mi are chosen to be 10^ GeV and 10^^ GeV respectively, and so T^ax ~ 2 x lO^^ Qgy. Ti k, Tmax- 
a is treated as constant and evaluated at Trch, with gi{Mz) obtained from aEM{Mz) = 1/128, 



s\n^ew{Mz) = 0.231, a^(Mz) = 0.119, and Mz = 91.2 GeV 



constrained than before. 

Comparison with numerical analysis: The total gravitino density, including that 
generated during reheating, has been obtained numerically in Ref. Our analytical 

derivation agrees well with their fit to the gravitino abundance. For a reheat temperature 
of 10^ GeV, both Eq. (F12) of Ref. 2^ and our Eq. ([28]) give Y^ = 2x lO^^^ indicating the 
robustness of our analysis. Furthermore, our analytical derivation allows one to appreciate 
various aspects of the gravitino abundance obtained in Ref. 29|. The fit to the gravitino 
abundance in Ref. j29|] gives an abundance dependent on Treh, but not also on Tmax as 
one might have expected. Our analysis above shows that the abundance generated during 
reheating does indeed depend on Tmax (see Eq. f l26l) ) but it also depends on the scale of 
inflation and Treh, and by manipulating the expressions relating these three quantities, as 
we have done, the dependence on Tmax cancels out. 



The dominant term in the fit of Ref. 



29| for the total gravitino abundance has the same 



functional form as that obtained by earlier analytical calculations that did not include grav- 
itino generation during reheating, namely, linear dependence on Treh- Our analysis indicates 
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FIG. 2: = Uq/s generated during the reheating era, the radiation dominated era, and the sum of 
contributions from both eras are shown as a function of the temperature T for the same parameters as 
in Fig. ([T]). Since Y^j in both eras is largely generated close to T^eh, a is evaluated at Tj-eh- The final 



value of y^^, is w 2 X lO^^^ 



that this is because gravitino generation during reheating also has a hnear dependence on 
Treh, just as in the post-reheating era. We emphasise that it would be improper to naively 
conclude that the hnear dependence is because the gravitino abundance generated during 
the radiation dominated era is dominant, since we have shown that the gravitino abundance 
generated in both eras differ only by a factor of 3. 

IV. CONCLUSION 



In conclusion, in this article we have calculated the gravitino abundance generated during 
reheating. We find that it is linearly proportional to the the reheat temperature Tj-eh, as in the 
standard calculation of gravitinos produced in the radiation dominated era after reheating. 
Further, we find that it is about 1/3 the number density of gravitinos generated in the 
radiation dominated era. Therefore this lowers the upper bound on Tj-eh from cosmological 
constraints on the gravitino number density by a factor of 4/3. This does not significantly 
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alter the viability of leptogenesis scenarios. 
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